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ABSTRACT 

We apply two different algorithms to search for mass segregation to a recent obser- 
vational census of the p Ophiuchi star forming region. Firstly, we apply the Amsr 
method, which compares the minimum spanning tree (MST) of a chosen subset of 
stars to MSTs of random subsets of stars in the cluster, and determine the mass seg- 
regation ratio, Amsr- Secondly, we apply the m — S method, which calculates the local 
stellar surface density around each star and determines the statistical significance of 
the average surface density for a chosen mass bin, compared to the average surface 
density in the whole cluster. Using both methods, we find no indication of mass seg- 
regation (normal or inverse) in the spatial distribution of stars and brown dwarfs in 
p Ophiuchi. Although p Ophiuchi suffers from high visual extinction, we show that 
a significant mass segregation signature would be detectable, albeit slightly diluted, 
despite dust obscuration of centrally located massive stars. 
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1 INTRODUCTION 

Most stars form in groups, clusters, and larger associations. 
In order to understand the star formation process, it is de- 
sirable to quantify the spatial distribution of stars in dif- 
ferent star forming regions, so that a clear picture of the 
formation and evolution of each region can be drawn. It is 
possible to measure the amount o f substructure in a region 
(e.g. by using the Q-parameter, ICartwright fc Whitworthl 
l2004h and to quantify the amount of ma ss segregation 
(e.g. the Amsr method, Allison et al. 2009a|, or the m — S 
method. [Maschberg er fc Clarke! 201 lj h Additionally, statis- 
tical methods can be applied to find clusters against a b ack- 
gro und field (e.g. iGutermuth et all 120091 : ISchmeial 1201 ll ) . 

[All ison et all l|2009al ) found that the amount of mass 
segregation in the ONC could be quantified by comparing 
the minimum spanning trees (MSTs) of chosen subsets of 
stars to the MSTs of random sets of stars. If the MST of the 
most massive stars is shorter than the MSTs of random sub- 
sets of cluster stars, then the clust er is mass segregated. The 
ONC is mass segregated (see also iHillenbrand fc Hartmannl 
Il998lh and the same signature was found by Sana et all 
j20ld i in Tru mpler 14. 

However. I Parker et alj (|201 lh found that the most mas- 
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sive stars in the Taurus association were 'inversely mass seg- 
regated', i.e. anti-clustered with respect to randomly cho- 
sen stars. Is mass seg regation therefore a dynamical pro- 
cess (as postulated by lAllison et al.ll2009bl h rather than a 
primordial outcome of star formation (in hydrodynamical 
simulations of star cluster formation, primordial mass seg- 
reg ation occurs as part of the competitive accretion proces s , 
e.g. iMaschberger et alJliuTol . iMasdiberger fc Clarke! 120111 )? 
To answer this question, we must first search for mass seg- 
regation in other young star forming regions, ideally using 
independent methods. 



In this paper we search for mass segregation in p Ophi- 
uchi. This cluster suffers heavily from differential extinc- 
tion, so an accurate and self-consistent determination of 
stellar masses is difficult. However, recent spectroscopic 
surveys jAlves de Oliveira et al.ll2012l : lErickson et aljfeoill ; 
ICeers et al.ll201ll : iMuzic et al.ll2012r ) have probed the low- 
mass end of the IMF and allowed a complete census of the 
cluster to be made. We describe the observational sample in 
Section [2l we describe the methods used to quantify mass 
segregation in Section O we present our results in Section U 
we discuss the results and the potential effects of extinction 
in Section [5] and we conclude in Section [6] 
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Figure 1. A map of p Ophiuchi showing the 255 objects in our datasct (restricted to the WIRCam field). The 20 least massive cluster 
members (masses up to 0.03 Mq) are shown by the blue crosses and the 20 most massive cluster members (masses down to 1.63 Mq) are 
shown by the large red dots. The solid lines indicate the extent of the WIRCam field. 



2 THE OBSERVATIONAL SAMPLE 

We adopted as a starting point for building the observa- 
tional sample the most re cent census of the p Ophiuch i 
core cluster as compiled in lAlves de Oliveira et al] l|2012l) . 
where a new population of brown dwarf members in the 
cluster was uncovered and analyzed with respect to the pre- 
viously known members. Their census includes only stars 
and brown dwarfs for which there is a spectral type classifi- 
cation and reliable membership confirmation. In short, the 
compilation takes all the spectroscopically co nfirmed mem- 
bers compiled for the cluster's review (see, IWilking et al.l 
I2008I . and refer ences therein) in the Handbook of Star 
Forming Regions (iReipurthl feoOS). adding to it the more re - 
cent spectrosco p ic res ul ts bvlAlves de O liveira ct al.l 1120101) 
iMcClure et ail (|2010t ): iGeers et all l|201lf ): Eric kson et al.l 
l|201ll ), totalling a list of 250 members where 208 have spec- 
tral types earlier than M6, and 42 have spectral types later 
or equal to M6. At the age of p Oph (~1 Myr), the evo- 
lutionary models of the Lyon group (|f3araffe et al.l 1 19981 ; 
IChabrieir et"aI]l200Ch w hen co mbined with the temperature 
scale of Luhman et al. J20031), plac e the substellar bound- 
ary at ~M6.25 ( Luhman et al.ll2007l ). We have added to this 
census two low mass stars a nd one brown dwarf recently con- 
firmed spectroscopically bv lMuzic et al.l (|2012l ). as well as 3 
members where the spectral type carries a larger error (but 



membership is confirmed) presented by Alves de Oliveira et 
al. (2012) but not included in their compilation. 

In the substellar regime, the spectroscopic follow-up o f 
the CFHT/ WIRCam survey (| Alves de Oliveira et al.ll2012t ) 
is nearly complete down to an extinction of 20 visual mag- 
nitudes (only 3 photometric candidates were not observed 
spectroscopically), and within the WIRCam mapped re- 
gion (interior to the solid lines in Fig. [TJ. In the stel- 
lar domain, a co nservative depth of 8 v isual magnitudes 
has been used bv lAlves de Oliveira et all l)2012l ) to define a 
complete sample with spectroscopic confirmation. To com- 
plete the database in the stellar domain at higher extinc- 
tions, we have included candidate members from X-ray sur- 
veys that still lack a spectroscopic confirmation. X-ray sur- 
veys of young stellar objects usually have low contamina- 
tion rates, in particular in clusters like p Oph where large 
amounts of extinction effectively block background sources. 
From the 51 X-ray sour ces which lack spectr oscopic confir- 
mation compiled in the IWilking et al. I (l2008h list of candi- 
date members (originally uncovered bv[lmanishi et al.ll200ll ; 
iGaene et al.ll2004l : lOzawa et al.ll2005l : IPillitteri et al.ll2010l ). 
43 are matched to a near-IR source, either in the 2MASS 
or the WIRCam catalogues. From the remaining 8 X-ray 
sources, 2 have an uncertain membership status (denoted 
as x? in IWilking et ail I200&D . and 6 are not detected in 
the J-band, and in any case are outside the photometric 
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completeness limits of the WIRCam surv ey (J=20.5 and 
#=18.9 mag. lAlves de Oliveira et alj|2012h . 

We have used the colour-colour diagram J— H vs. H — K 
to dere dden each of the X-ray sources along the extinction 
vector (|Rieke fc Lebofskvl 19851 ) and estimate a spectral type 
by comparing their near-IR photosph eric colours to those 
characteristic of young stellar objects l|Luhman et aljfeoid . 
their Table 13). This method could not be applied for 9 
sources which have strong IR excess and therefore their po- 
sition on the colour-colour diagram is likely to be affected 
by the contribution of the disk, and to 4 sources which al- 
though classified as Class 111 fall in a region of the diagram 
where the colours of young stellar objects (~M9 to early L) 
increase nearly parallel to the reddening vector, and there- 
fore any solution is degenerate. We estimated spectral types 
(~K4 to ~M8) for 30 sources using this method. 

To estimate the masses of the members of the cluster, 
we first convert spectral types to temperature , adopting the 
temperature scale from ISchmidt-Kalerl (ll982T) for stars ear- 
lier than MO, and the scale from Luhman et al.l l|2003l ) for 
sources with spectral type between MO and M9.5. For the 
L dwa rfs, we applied the scale proposed by lLodieu et al.l 
(2008) extrapolated to the L4 spectral type. Masses were 
deriv e d from the 1 Myr ev o lutionary models (|Baraffe et al.l 
1 19981 ; IChabrier et~all 120001 ; ISiess et all 120001) according to 
each target's effective temperature. 

Because we derive our masses from temperatures, many 
stars are assigned the same mass from the stellar models. 
This is potentially a problem for our mass segregation algo- 
rithm, as 25 stars at the peak of the IMF may be assigned 
the same mass. We therefore apply a small amount of ran- 
dom noise to each mass, thereby making each value unique. 

In Fig.[T]we show the 255 objects in our dataset (which 
we have restricted to the WIRCam field). We plot the 20 
most massive stars (masses down to 1.63 Mq) as the large 
red points, and the 20 least massive objects (masses up to 
0.03 Mq) as the blue crosses. 



2.1 Spatially and extinction-limited sample 

A major caveat in studying a representative sample of 
the p Ophiuchi population is the variable extinction across 
the cluster. To attempt to correct our methods for this, 
we also examined a spatially and extinction limited sam- 
ple of objects. We selected from the original data base 
all sources that had in the colour-magnitude diagram 
H vs. J—H an Ay<20 mag (see, for example, Fig. 9 in 
lAlves de Oliveira et al.ll2012h . The data base limited both 
spatially and to an extinction of 20 visual magnitudes con- 
tains 205 members. Though the masses of the X-ray mem- 
bers determined from photometry are likely to carry a large 
uncertainty, they represent only ~11% of this sample, and 
should nevertheless reflect in relative terms the relation be- 
tween the real masses. 



3 METHOD 

In this section we outline the two methods we use to look for 
mass segregation signatures in the dat a, namely the Amsr 
ratio pioneered bv lAllison et aL (12009a ) and the m— £ distr i- 
bution, recently proposed by [Maschberger fc Clarkel |201ll ). 



3.1 The Amsr mass segregation ratio 

We first quantify any mass segregatio n present in the cluste r 
by using the Amsr ratio introduced bv lAllison et al.l (l2009al ). 
This constructs a minimum spanning tree (MST) between a 
chosen subset of stars and then compares this MST to the 
average MST length of many random subsets. 

The MST of a set of points is the path connect- 
ing all the points via the shortest possible path l ength 
but which contains no close d loops (e.g. iPriml 1 19571 : 
ICartwright fc Whitwo"rthll2004t). 

We use the algorithm of lPriml (|l957t l to construct MSTs 
in our dataset. We first make an ordered list of the separa- 
tions between all possible pairs of starfl Stars are then con- 
nected together in 'nodes', starting with the shortest separa- 
tions and proceeding through the list in order of increasing 
separation, forming new nodes if the formation of the node 
does not result in a closed loop. 

We find the MST of the Nmst stars in the chosen sub- 
set and compare this to the MST of sets of A?mst random 
stars in the cluster. If the length of the MST of the chosen 
subset is shorter than the average length of the MSTs for 
the random stars then the subset has a more concentrated 
distribution and is said to be mass segregated. Conversely, 
if the MST length of the chosen subset is longer than the 
average MST length, then the subset has a less concentrated 
dist ribution, and is sai d to be inversely mass segregated (see 
e.g. iParker et al1l201ll ). Alternatively, if the MST length of 
the chosen subset is equal to the random MST length, we 
can conclude that no mass segregation is present. 

By taking the ratio of the average (mean) random MST 
length to the subset MST length, a quantitative measure of 
the degree of mass segregation (normal or inverse) can be 
obtained. We first determine the subset MST length, Z au b se t. 
We then determine the average length of sets of A'mst ran- 
dom stars each time, (/average). There is a dispersion asso- 
ciated with the average length of random MSTs, which is 
roughly Gaussian and can be quantified as the standard de- 
viation of the lengths (/average) ± ""average • However, we con- 
servatively estimate the lower (upper) uncertainty as the 
MST length which lies 1/6 (5/6) of the way through an or- 
dered list of all the random lengths (corresponding to a 66 
per cent deviation from the median value, (/average))- This 
determination prevents a single outlying object from heav- 
ily influencing the uncertainty. We can now define the 'mass 
segregation ratio' (Amsr) as the ratio between the average 
random MST pathlength and that of a chosen subset, or 
mass range of objects: 



At, 



(1 



average 



^ + (J 5/6/ Z subsct 
1/6/^subset 



(1) 



A Amsr of ~ 1 shows that the stars in the chosen subset are 
distributed in the same way as all the other stars, whereas 
Amsr > 1 indicates mass segregation and Amsr < 1 indi- 
cates inverse mass segregation, i.e. the chosen subset is more 
sparsely distribut ed than the other star s. 

As noted by lAllison et all (|2009al ). the MST method 
gives a quantitative measure of mass segregation with an 

1 From this point onwards, when referring in general to 'stars' in 
the cluster, we mean 'stars and brown dwarfs', as we are including 
all the objects in the observational sample. 
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associated significance and it does not rely on denning the 
centre of a cluster. 

There are several s u btle variations of Amsr- 
lOlczak. Spurzem fc Hennind i|201ll ) propose using the 
geometric mean to red u ce the spread in uncertainties, and 
IMaschberger fc Clarke! (|201lT ) propose using the median 
MST length to reduce the effects of outliers from influencing 
the results. However, in the subsequent analysis we will 
adopt the original Amsr from Allison. 



3.2 The m - E distribution 



1 .63 



.27 0.66 0.38 0.20 0.12 0.03 



Recently, IMaschberger fc Clarke! (|201l! ) proposed a method 
to analyse mass segregation which measures the distribution 
of local stellar surface density, E, as a function of stellar 
mass. We calculate t he local stellar sur f ace d ensity follow- 
ing the prescription of lCasertano fc Hud (|l985l ). modified to 
account for the analysis in projection. For an individual star 
the local stellar surface density is given by 



N ■ 



(2) 



where rjv is the distance to the N th nearest neighbouring 
star (we adopt N = 10 throughout this work). 

If there is mass segregation, massive stars are concen- 
trated in the central, dense region of a cluster and thus 
should have higher values of E. This can be seen in a plot of 
E versus mass, showing all stars and highlighting outliers. 
Trends in the m — E plot can be shown by the moving 
average (or median) of a subset, E su b sc t, compared to the 
average (median) of the whole sample, E a ii. The signature 
of mass segregation is then E su b S ot > E a ii, and of inverse 
mass segregation E su b sc t < E a n. The statistical significance 
of mass segregation can be established with a two-sample 
Kolmogorov-Smirnov test of the E values of the subset 
against the E values of the rest. 

Note that there are many more ways of defining mass 
segregation. For instance, one can choose a cluster centre 
and m easure the mass function as a function o f radial dis- 
tance (|Gouliermis et al.ll200j ; ISabbi et al.ll2008l ). or the dis- 
tance of the most massive star(s) from the cluster centre 
compared to t he average distance of low-mass stars to the 
cluster centre (|Kirk fc Myers! l201ll ). Both methods rely on 
determining the centre of the cluster or association, which 
in the case of low-number clusters with substructure is non- 
trivial and is virtually impossible in the case of a hig hly 
substructured region such as Taurus (I Parker et al.|[201lf ). 



4 RESULTS 

In this section we present the results of our Amsr analysis, 
followed by the m — E distribution. We then discuss the 
effects of extinction on the results. 



4.1 A„ 



for high mass stars 



In Fig. [2] we show the evolution of Amsr as a function of 
the number of stars in an MST, ./Vmst for the most massive 
stars in the cluster. We increase the number of stars in the 
MST in steps of 6, which is a compromise between a high 




tftfTfififinnnnt 



Figure 2. The evolution of the mass segregation ratio, AmsRi 
for the A^mst most massive stars in our dataset and we also indi- 
cate the lowest mass star, rriL within the NyisT- Error bars show 
the 1/6 and 5/6 percentile values from the median, as described 
in the text. The dashed line indicates Amsr = lj i- e - 110 mass 
segregation. 



enough resolution to pick out structure between different 
mass regimes, and a low enough resolution so that we do not 
add noise to the plot. The first subset compares the MST of 
the 20 most massive stars to the median of many different 
random sets of 20 stars, and the second subset is the 26 
most massive stars compared to the median of random sets 
of 26 stars, and so on. On the top axis we also indicate the 
mass of the least massive star within that value of ./Vmst , at 
regular intervals. 

Fig. [2] we see that there is no clear mass segregation 
signature (normal or inverse) in the most massive stars in 
the cluster (the most massive 20 stars are indicated by 
the large red points in Fig. [TJ). The 20 most massive stars 
(with masses above 1.63 Mq) have a mass segregation ratio 
Amsr = 0.891q'?3, which does deviate from Amsr = 1 (indi- 
cating slight inverse mass segregation), but because the 26 
most massive stars are consistent with Amsr = 1, this result 
is not particularly significant. 



4.2 Amsr for low mass stars 

In Fig.|3]we show the evolution of Amsr as a function of the 
number of stars in an MST, ./Vmst for the least massive stars 
in the cluster. We begin by constructing an MST with the 
20 least massive objects in the cluster, and then increasing 
the number of objects in the MST by 6 at each stage. On 
the top axis we now indicate the mass of the most massive 
star within the ./Vmst subset. 

We see that the least massive objects do not show any 
strong mass segregation signature, and (within the uncer- 
tainties) are consistent with Amsr = 1. 



4.3 The m — E distribution 

We show the m — E distribution for the stars in our dataset 
in Fig. [4] The upper (black) dashed line is the mean E value 
for the whole cluster, and the lower (blue) dashed line is the 
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Figure 3. The evolution of the mass segregation ratio, Amsr> f° r 
the A^mst least massive stars in our dataset and we also indicate 
the highest mass star, mji within the Nmst- Error bars show 
the 1/6 and 5/6 percentile values from the median, as described 
in the text. The dashed line indicates Amsr = li i- e - no mass 
segregation. 



Figure 4. The m — S distribution for the stars in our dataset. We 
plot the local surface density for each star against its mass. We 
then determine the median (blue lines) and mean (black lines) 
S for the entire cluster (the dashed lines) and for the 50 least 
massive, and 50 most massive stars in the cluster (the solid lines). 



median value. We also show the mean and median E values 
for the 50 most massive stars (on the righthand side) and 
the 50 least massive stars (on the lefthand side) by the solid 
lines. 

The plot shows that the mean and median E values of 
the lowest mass objects in the cluster are marginally higher 
than for the whole sample. The p-values of a two-sample 
KS test (E of low-mass stars versus the entire cluster) are 
p = 0.21 (20 least massive) and p = 0.51 (50 least massive). 
Usually, these would need to be smaller than p — 0.05 at a 
significance level corresponding to 2a, in order to reject the 
hypothesis of "no mass segregation" . Thus, the lowest mass 
objects are not mass segregated. 

The most massive stars lie at slightly lower E values 
compared to the whole cluster, suggesting inverse mass seg- 
regation. Here the p- values are p — 0.17 and p = 0.70 for 
the 50 and 20 most massive stars, respectively. Again, this 
does not indicate any significant deviation of the spatial dis- 
tribution of the massive stars from the spatial distribution 
of the other stars. The 50 most massive stars are inversely 
mass segregated, similar to the Amsr. results, but only at a 
la level. This is not the case for the 20 most massive stars, 
where no inverse mass segregation can be concluded. Given 
the small n and the rather weak signature for Amsr this 
result can be deemed compatible with the Amsr- 



4.4 Extinction-limited sample 

A major caveat in determining the spatial distribution of a 
sample of objects in p Ophiuchi is the variable extinction 
across the cluster. As a check that our results do not change 
when an extinction limit is imposed on the data, we apply 
an A v limit of 20 mag and then repeat the MST and m — 
E analysis on this extinction-limited sample. We find no 
discernible difference to the results in either case; i.e. there 
is no clear mass segregation signature in either the high- or 
low-mass objects in the cluster. 



5 DISCUSSION 

The results presented in Section [4] show that there is no ev- 
idence of mass segregation in that the most massive stars 
are not centrally con centrated, as they a re in e.g. the ONC 
jAllison et alj l2009al ) and Trumpler 14 (|Sana et alj|2010h . 
This could indicate that mass segregation may be a dynam- 
ical process, rather than a primordial outcome of star for- 
mation, but a study of more star forming regions is required 
to substantiate this hypothesis. 

In this dynamical scenario, the massive stars form at 
random locations in a substructured cluster, and then a sub- 
virial collapse facili tates mass segregati on on a very short 
timescale ( ~ 1 Myr. lAllison et al.l 2009bl). o Oph is not sub- 
structured l|Cartwright fc Whitworth |2004| ). but may have 
been at earlier ages. If it was substructured at earlier ages, 
this has not facilitated dynamical mass segregation in this 
cluster. 



5.1 Extinction 

The high level of extinction makes observing objects in 
p Oph challenging, and it is possible that even with our 
extinction-limited sample, some stars are still hidden in the 
centre of the cluster. In such a scenario, unobserved high 
mass stars could reside in the central regions, and any mass 
segregation of such stars would not be observed. In this case, 
both our mass segregation-finding algorithms would erro- 
neously give a null-result, similar to those described in the 
previous Section. 

Here, we conduct a simple numerical experiment to de- 
termine how much a mass segregation signature could be 
diluted by high levels of extinction, such as that present 
in p Op h. We distribut e 360 stars randomly in a Plummer 
sphere |Plunimerlll91ll ). with a half-number radius of 1 pc 
accord ing to the prescription in lAarseth. Henon fc Wielenl 
( 1974), and as sign masses (again at random) from a 3-part 
lKroui.nl (|2002h IMF of the form: 
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(a) No mass segregation 




(d) Mass segregation 




2.929 0.526 0.290 0.158 0.095 0.064 0.034 



OffifJllpftlfifi* 



(b) No mass segregation 



m L (M s ) 

2.929 0.526 0.290 0.158 0.095 0.064 0.034 



i(» if > 



(e) Mass segregation 

m L (M Q ) 

2.282 0.467 0.290 0.158 0.086 0.053 0.027 



! 



(h) Mass segregation with extinction 




Moss (y e ) 
(c) No mass segregation 



Moss (M G ) 

(f) Mass segregation 



Mass (M G ) 



(i) Mass segregation with extinction 



(g) Mass segregation with extinction 

Figure 5. Results for a fake cluster without mass segregation (panels a — c), with mass segregation (d - f), with mass segregation but 
some stars removed due to high extinction (g - i) . In each case we show the morphology of the cluster (the most massive stars are shown 
by the large red dots), the plot of Amsr f° r the high mass stars, and the m — S plot. The extent of the extinction cloud is shown by the 
black circle, in the m — X plots the median £ value for the most massive stars is shown by the solid red line, and the median £ value for 
the whole cluster is shown by the dashed line. 



N(M) 




m < M/M Q ^ mi . 
mi < M/Mq ^ m 2 , 
m 2 < M/Mq ^ m 3 . 



(3) 



and we choose mo = 0.02 Mq, mi = 0.1 Mq, m 2 = 0.5 Mq, 
and m-3 = 10 Mq. In Fig. 5(a) we show the morphology 



of this cluster, with the 40 most massive stars shown by 
the red points. If we determine Amsr for the most massive 
stars (in steps of 6 objects) we see that this cluster is not 
mass segregated, with Amsr — 1 throughout (Fig. 5(b) I. 
The m — E algorithm also shows no significant differences 
between the 40 most massive stars and the cluster as a whole 
(the solid red line and the dashed line, respectively, shown 
in Fig. 5(c) I. 



We apply a simple mass segregation algorithm to the 
Plummer sphere by swapping the positions of the 40 most 
central stars with the positions of the 40 most massive stars 
(we choose 40 stars to clearly demonstrate the effects of 
extinction in Fig. [5] but the results are similar for the 20 
most massive stars). We show the new spatial configuration 



of the massive stars in the cluster in Fig. 5(d) Several of 



the most massive stars are originally within the sample of 
the 40 most central stars, and end up (randomly) being as- 
signed positions outside the central core. In one sense, such 
a configuration is perhaps more realistic than if the 40 most 
massive stars were also the 40 most central; in a real clus- 
ter dynamical interactions between the central stars would 
likely eject one or two of the massive stars. 
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In Fig. 5(e) we show the evolution of Amsr as a func- 
tion of the number of stars in the MST. The effect of ar- 
tificially mass segregating the cluster is clearly seen, with 
Amsr = 18.6lg°7 2 . The cluster shows significant mass seg- 
regation down to the 40 th most massive star, which has 
Amsr = 3.6lg|. Similarly, the m — E method also shows 
that the cluster is mass segregated; in Fig. 5(f) we show 



the median surface density of the entire cluster by the 
dashed line (E = 6.35 stars pc -2 ) and the median surface 
density of the 40 most massive stars by the red solid line 
(E = 10.76 stars pc -2 ). A two-sample KS test returns a p- 
value of < 10 -8 that the two distributions could be drawn 
from the same parent population. 

We now assign a power-law extinction to the fake clus- 
ter, from the centre out to a radius of 5 pc (denoted by the 
circle in Figs 5(a) 5(d) and 5(g) I. We then assign an A v 
value to each star using the following formula: 



A v (r) = 20 



r 
5pc 



5/3' 



(4) 



where r is the position of the star with respect to the cluster 
centre and \r\ is the modulus of its vector. To account for 
projection effects along the line of sight, we double A v (r) if 
the z-component of the vector r is negative. Therefore, in 
the central regions of the cluster, the A v value can range 
between ~ 10 — 40. We then remove all stars with A v > 20, 
leaving a total of 193 stars. In Fig. |5(gj] we show the spatial 
distribution of the remaining objects. 

Once again, we calculate Amsr for the remaining ob- 
jects, and Fig. 



5(h) 



shows that the mass segregation sig- 
nature is still observable, although to a lesser extent due 
to the removal of several of the most massive stars in the 
cluster. The peak value is now Amsr = 6.4^3'g, but the 
plot still shows the same morphology as the non-extinction- 
limited data-sample. Furthermore, the m — E method also 



shows that the cluster is still mass segregated; in Fig. 5(i) 



we show the median surface density of the entire cluster by 
the dashed line (E = 4.50 stars pc -2 ) and the median sur- 
face density of the 40 most massive stars by the red solid 
line (E = 7.34 stars pc~ 2 ). A two-sample KS test returns a 
p-value of < 10~ 2 that the two distributions could be drawn 
from the same parent population. 

We have demonstrated with a simple model for extinc- 
tion that the two mass segregation finding algorithms could 
still determine whether a cluster suffering from extinction 
is significantly mass segregated or not. The results suggest 
that the actual dataset, whilst possibly lacking some clus- 
ter members due to obscuration, is likely reflecting the true 
spatial distribution of stars and brown dwarfs in p Oph. 



6 CONCLUSIONS 

We have used an observational census of p Ophiuchi, which 
was recently enhanced by several surveys probing the sub- 
stellar domain of the IMF, to search for possible mass seg- 
regation signatures in the spatial distribution of stars and 
brown dwarfs in this cluster. 

We have utilised two different algorithm s. Firstly, we 
used the Amsr technique (| Allison et alj|2009af ). which com- 
pares the minimum spanning tree (MST) of a chosen subset 



of stars, to the MSTs of randomly chosen stars in the clus- 
ter. If the MST length of a chosen subset is shorter than 
the MST length of the random objects, then the cluster is 
mass segregated. Secondly, we have used the m — E plot, 
which compares the local surface density surrounding mas- 
sive stars to the the average surface density of all of the stars 
in the cluster. By this definition, a cluster is mass segregated 
if the massive stars have a significantly higher than average 
surface density. Our conclusions are as follows: 

(i) The Amsr technique finds that the most massive 
stars show hints of being inversely mass segregated, with 
Amsr = 0.89±g;?| for the 20 most massive stars. However, 
Amsr is consistent with there being no mass segregation 
of the 26 most massive stars, and so on. The least massive 
stars show no clear deviation from Amsr = 1- 

(ii) The m — E distribution also suggests that the most 
massive stars may be inversely mass segregated (but with 
no strong statistical significance), and with no difference in 
the distribution of low-mass stars compared to the cluster 
average. 

(iii) The high levels of extinction in p Oph may mean that 
some members are missing from the dataset. However, 
we have demonstrated that a significant difference in the 
spatial distribution of a group of objects would still be 
found by both the Amsr and m — E methods. 

In order to understand the star formation process in 
different clusters, we suggest applying both mass segregation 
algorithms in tandem to build up a census of the spatial 
distribution of stars in different star forming regions. 
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